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PROGRAMMED INSTRUCTION - PROJECT 1000 


W. J. Lnenicka 
Georgia Institute of Technology 


INTRODUCTION OF PROGRAMMED INSTRUCTION 


Ideas, methods and techniques which are new or different 
in education are apt to create controversy and to be 
resisted. Programmed instruction is no exception. 


Debate and argument will not solve the questions of 
whether to use, or how to use, programmed instruction. 
The only sure way to find out whether or not programmed 
instruction can be beneficial is to try out, improve 
and evaluate programs in a typical educational setting. 


Administrators should assume the responsibility for the 
introduction and evaluation of programmed instruction, 
just as they do for all other kinds of instructional 
innovations. A cautious approach is important when any 
school considers committing its curriculum, or any part 
of it, to programmed learning. Careful attention must 
be paid to preparing the teaching staff, preparing the 
students, assessing the costs, providing administrative 
support, and selecting and preparing programs. 


The preferred way to introduce programming into a school 
is to encourage an interested, informed, capable teacher 
to try a program on a limited basis. Programs covering 
only one topic can be tried without disrupting current 
school practices. For the informed teacher with some 
experience, programmed material covering an entire 
course or term can be used. Of course, it should be 
obvious that programmed material must be well organized 
and well written to be useful. Teachers who would not 
use a second-rate textbook should not feel that poor 
quality programmed material is suitable for student 
usage. 


All programs designed for classroom use need to be 
supplemented, They are not materials designed to 
replace teachers, to avoid the need for discussions, or 
to replace all contact between students and professors. 
The lecture is not dead, discussion is not dead, class- 
room contact is not dead, and those few who make such 
claims invariably do so by lecturing to their audience 
without seeming to recognize the contradiction. 


PREPARING TEACHERS TO USE PROGRAMMED INSTRUCTION 


The reaction of teachers to programmed instruction is 
varied. Some regard programming as useless, or as a 
means to water down the rigorous courses in the curricu- 
lum. Others view the introduction of programmed material 
as a device to cut down on their student contact hours 
(reduce their teaching load); still others respond with 
unbounded enthusiasm and view programming as a method 
that will revolutionize education, Most knowledgeable 
instructors experienced in programmed instruction 
operate somewhere between the extremes. Some teachers 
are favorably impressed but others respond negatively 

to programming. Why? Perhaps because programs are new, 
and published programmed texts on engineering subjects 
are scarce. Teachers are untrained in the use of pro- 
grams, schools are not geared to their use, particularly 
when they are used in self-paced ways, and some faculty 
members are openly hostile to departure from traditional 
practice. In any event, there are two key factors in 
the successful use of programmed instruction. One key 
is the programmed material; the other key is the teacher. 


Provision should be made for preparing the teaching staff 
to use programmed materials. The minimum training that 
teachers need must produce an understanding to the basic 
learning theory relative to programming, the role that 
specific objectives play in program development, the 
sequencing of content, the basic principles of frame 
writing, the uses to be made of student responses, the 
interpretation of test data, and the correlation of 
programs with other instructional activities. 


Teachers should continue to be the directors and the 
skilled guides to the learning process. Programs are 
no threat to a teacher's existence, Instead, they are 
a tool which, if properly prepared and utilized, can 
help students reach selected educational objectives 
efficiently and effectively. Teachers and students will 
enjoy using programs which challenge their competence 
and, at the same time, aid in effective learning. 


Teachers who will be required to use programs should 
have the opportunity to develop a basic understanding 
of programmed learning. The opportunity can be provided 
by enrollment in courses, seminars, workshops, or inten- 
sive study of selected written materials. It is a mis- 
take to assign teachers with no prior knowledge or 
experience with programmed methods to teach programmed 
courses. 


Teachers should be capable of taking a body of program- 
med material and selecting from it programs for 

specific purposes, Most teachers do not hesitate to 
select appropriate material in textbooks for student 
assignment; yet many who are faced with programmed 
material prepared by others seem to think that all of 

it must be assigned to the students "because it is 
there." Attitudes such as this reflect not only inexper- 
ience but also relate to a basic lack of knowledge on 
the part of the teacher. 


Some students react negatively to programs and do poorly 
with them. It may be disinterest, a reading handicap, 
inability to study alone, or for reasons not understood, 
Usually a diversity of programmed material sufficient 

to appeal to these students is not available; teachers 
should advise them to enroll in traditional classes, if 
possible. 


Some of the better students become bored, unless they 
are permitted to pace themselves in their progress 
through the course. When self-pacing is permitted, the 
superior students regard finishing the course early as 

a challenge, and usually take off running. Most are 
extremely proud and satisfied if they can finish several 
weeks early. Naturally, a course does not need to be 
programmed to be self-paced, but, on the other hand, the 
very nature of programmed material lends itself to self- 
pacing. 


Since experience has clearly shown that a group of 
learners will vary widely in their speed of learning 

and progress through the course, the teacher will soon 
find that the students are spread through widely varying 
subject matter areas. This means, as a practical matter, 
that the teacher must know the whole subject thoroughly 
before the course begins, This feature of programmed 
instruction will frighten teachers who try to stay only 
one day ahead of their students. 


In addition, the teacher will need a considerable depth 
of understanding in the subject matter. Because the 
better students go faster, they will need more subject 
matter, harder problems, or higher level subjects, 
particularly if self-pacing through the course is not 
permitted. Teachers with deficiencies in subject matter, 
or who lack confidence in their ability to handle new 


and different situations, should not attempt programmed 
instruction, 


Teachers need to be quick to recognize those students 
who are not progressing satisfactorily, and to attempt 
to find a solution to the problem. As in any kind of 
teaching situation, it is easier to detect the existence 
of a problem than it is to do something positive about 
it. Students may complain of boredom, lack of under- 
standing of the frames, and of excessive slowness. They 
may be unable to concentrate and they may fail the tests. 
The failure to learn may be due to causes within the 
learner, the program or the teacher. Programmed learn- 
ing needs the support of continued and frequent testing, 
with aids and procedures for remedial learning. 


PROGRAMS AS HOMEWORK 


Programmed material designed for use in a university 
setting is designed primarily for students' use at a 
time and place chosen by the student. In other words, 
programs are designed for use as homework. Failure to 
recognize this fundamental fact has led some to predict 
that teachers will become learning resource managers, 
whatever those are. Teachers in engineering who use 
programmed material will need to continue doing what 
they have always done--clarify concepts, provide motiva- 
tion, monitor student progress and evaluate the results. 


The use of programmed material, properly prepared and 
utilized, increases the student's efficiency in doing 
his homework, and makes it possible for him to get more 
information on his own. That is all! 


If a teacher recognizes when and where programmed 
material is being used and for what it is being used, 

he can view the situation in its proper perspective. 
Traditionally, homework assignments have frustrated 
learners, encouraged copying, developed in students a 
dislike for school, and created antagonism among students, 
parents and techers. The careful construction of 
programmed frames, with the immediate reinforcement pro- 
vided for the successful accomplishment of small tasks, 
lends itself to the development of better attitudes 
about homework, besides increasing the efficiency of 

the process. 


Good homework should be individualized, require indepen- 
dent study, relate directly to the course objectives, 
challenge the understanding and skill already possessed 
by the learner, and motivate the students to learn, 


In programming, we find a method which allows individual 
progress at a rate controlled by the learner. We find 
that independent study is required by the nature of the 
material; it would be impractical to have someone else 
respond to the frames. In addition, the programs are 
designed to produce successful responses rather than 
unsuccessful ones, removing some of the motivation to 
cheat. If self-pacing in the class is permitted, the 
varying abilities of the students soon result in a 
considerable scatter over the length of the course, 
making independent study the practical way to get the 
work done. 


Programs by definition are carefully designed to produce 
terminal behaviors that are clearly associated with 
course objectives. Through careful construction of 
programs the material relates directly to course 
objectives, thereby eliminating irrelevant material and 
unnecessary duplication. 


In programming, we find material that should challenge 
the understandings and skills already possessed by the 
learners. The students should be required to think, to 
work up to their abilities, and to use as many of the 
skills they possess as is possible or reasonable. Pro- 
grams, because of their careful sequencing and active 
Tesponses, keep the learner steadily progressing and 
achieving. The student is continually being challenged 
by new material. The danger of the learner getting into 
problems beyond his grasp is reduced because of the 
small steps, small information increments, and the 
logical sequencing of the material. 


Programmed material can motivate the learner favorably. 
Homework must compete with television, movies, games, 
sleep and other distractions, Homework assignments 
should be stimulating and worthwhile. Programs help 
to meet these requirements, and provide success which 
encourages progress. 


The use of programmed material permits the student a 
more effective means of acquiring new basic information 
and understanding. The teacher has more class time to 
clarify concepts, to provide more depth, and to increase 
the effectiveness of the class in problem-solving 
techniques. 


Many administrators and teachers are concerned with the 
extent to which a curriculum should be programmed. The 
answer is to use programs whenever their effectiveness 
can be demonstrated. The first step in finding out 
whether or not a program will work is to try one. 


PROGRAMMED INSTRUCTION IN MECHANICS AT GEORGIA TECH 


During the past few years the innovations in teaching 
methods employed in engineering mechanics courses at 
Georgia Tech have consisted of the use of programmed 
textual material in self-study courses. The particular 
courses that are currently utilizing programmed material 
and being conducted on a self-study basis are the under- 
graduate courses in statics and dynamics, 


Dynamics 


The possibility of teaching mechanics courses by use of 
programmed instruction material had been seriously 
considered for several years, The first attempt was 
made in the area of undergraduate dynamics. Since the 
summer of 1970, one or two sections of dynamics have 
been taught every quarter from a programmed text. 

The text used was A Programmed Introduction to Dynamics 
by Clyde E, Work; it was commercially available from 
McGraw-Hill Book Company. Other sections of the same 
course are still taught using a conventional text and 
traditional teaching methods. In this way it has been 
possible to compare the performance of students in the 
regular courses. 


The instructor has been granted freedom by the adminis- 
tration to experiment with various ways of conducting 
the course. As a result, the operation of the course 
has proceeded from one extreme to the other. At the 
one extreme there are no formal class sessions; the 
students come to the instructor's office for individual 
assistance when needed. At the other extreme there is 
a fairly rigid schedule of lectures and help sessions. 
Experiments with no formal lectures at all have not 
proved satisfactory, Many students fail to grasp the 
concepts readily on their own, but yet do not seek the 
instructor's help probably because they fail to realize 
they need assistance until after the exams. On the 
other hand, the schedule of frequent lectures with the 
programmed text serving much the same function as a 
nonprogrammed text fails to take full advantage of the 
programmed material, in that programmed material should 
allow the student to understand the concepts and 
problems with less assistance than is normally required 
in a course utilizing a conventional text. As a result 
of experimentation with various systems it has been 
found that for the five-hour, one-quarter course in 
dynamics the instructor needs to meet the class formally 
twice each week in order to clarify concepts, to discuss 


and illustrate problem-solving techniques, and to pro- 
vide motivation, 


The student response to this innovation in teaching 
dynamics has been favorable as is evidenced by what they 
tell the faculty, by their performance in the course, 
any by what they tell their friends, The five-hour, one 
quarter course in dynamics is no longer required in any 
engineering curriculum, but many students are continuing 
to take dynamics in this fashion either in place of the 
current three-hour, one quarter required course, or as 
an elective. For example, in the spring quarter of 

1973 there were 62 students enrolled in the course, even 
though their curricula required only a three-hour 
dynamics course for credit in this area of study. The 
programmed sections regularly reach capacity early in 
the registration period. 


Student performance on examinations has been similar to 
that of other sections conducted using nonprogrammed 
texts with a regular schedule of five lectures per week, 
A typical class averages about 2.4 on a 4.0 grade scale 
either way. Students have claimed that the programmed 
course requires less work than the regular course. It 
is expected that this would be the case since, in the 
regular course, much of the student's time is wasted in 
the instructor's attempt to simultaneously satisfy the 
diverse needs of a variety of students. In the pro- 
grammed course the individual can easily adapt the 
course to fill his individual needs. 


Because of the obvious saving in instructor time, more 
students can be handled; this also pleases the adminis- 
tration. The student-credit-hour load of the instructor 
teaching programmed dynamics is consistently among the 
higher ones in the department (about 25% higher than 
before). 


Statics 


Until quite recently no programmed text for statics 
comparable to the programmed dynamics text cited earlier 
was available from a commercial publisher. Thus, the 
initiation of programmed instruction in statics was 
necessarily delayed. However, favorable experience with 
the dynamics course led to the decision early in 1971 
to develop Georgia Tech's own statics text using a 
programmed format, Three faculty members financed by 
Institute and NSF funds devoted the equivalent of about 
one man year to the development of suitable material. 
The result was a programmed statics text in two volumes, 


printed on campus, for use beginning with the fall 
quarter of 1971. 


Five sections with a total of 142 students were taught 
from the programmed material during the 1971 fall quarter; 
more than 350 other students were enrolled in statics 
sections which were being taught in the traditional way, 
using Vector Mechanics for Engineers by Beer and Johnston 
as the text. The faculty arise teaching the programmed 
sections were free to administer the sections as they 
chose and were encouraged to experiment with innovations 
of any reasonable kind in the administration of the 
course. It was decided to divide the material in four 
main parts and to have an exam on each of these parts. 
The division was as follows: 
1. Moments, forces, free-body diagrams, equivalent 
force systems. 
2. Centroids and centers of mass, equilibrium 
concepts. 
3. Beams, trusses, frames and machines, 
4. Equilibrium in three dimensions, friction. 


Seven dates for giving examinations to students in all 
sections as a group were chosen. Any student who wanted 
to repeat an exam was permitted to do so, within the 
limits imposed by the seven dates. 


As a result of this effort, student performance was some- 
what higher in the programmed sections than in those 
using a conventional method. This was due to the self- 
pacing features of the programmed instruction sections, 
which allowed time to repeat examinations in a limited 
number of cases until the material was mastered. 


During the winter quarter of 1972 five programmed sections 
were offered under four different instructors. During 
registration, a fifth section was added to the four 
originally scheduled but, even so, several dozen students 
who asked to enroll in the programmed statics were turned 
down due to lack of space. The final figures showed 169 
in the programmed sections and 301 in the sections which 
used traditional methods. A group final examination was 
scheduled at the close of the winter quarter for all 
statics sections. The performance of all sections was 
similar. For the spring quarter, enrollment in the 
programmed sections was cut off when the balance of the 
500 texts which had been printed for use in 1971-72 was 
committed at registration. 


A few of the observations made during the conduct of the 
programmed sections in 1971 will now be discussed. First, 
as might be expected, there was a general tendency on 

the part of the student when required to budget his own 
time to get behind in the course and then try to do too 
much during the last days of the quarter, Also, without 
the urgency of day-by-day assignments, students tended 


to study statics only when the demands of their other 
courses had been met, and in many cases this amounted 
to a few days of concentrated study immediately prior 
to taking an examination. Secondly, with the course 
divided into four parts, a student could deceive himself 
into thinking that he understood one of these parts 
quite well and then find out otherwise as a result of 
the examination. In many cases, by that time the 
student was sufficiently behind that it was very diffi- 
cult to take action to remedy the situation, since he 
had to not only go back and repeat the material he had 
failed but at the same time proceed forward into the 
new material. 


Based on the experience with using the text during the 
1971 school year, it was decided during the summer of 
1972 to revise the text rather substantially, clari- 
fying some of the more difficult concepts, adding more 
problems where needed, correcting errors in the original 
text, and making other changes to more properly adapt 
the text to fulfill its objectives. 


Prior to the fall quarter of 1972 several changes were 
made in the administration of the programmed sections. 
The course was divided into ten parts, each part 
corresponding to a unit in the text. The division was 
as follows: 

1. Vectors 

2. Forces and moments 


3. Equivalent force systems 

4. Centroids and centers of mass 
S$. Equilibrium concepts 

6, Frames and machines 

7. Trusses 

8. Beams 

9. Equilibrium in three dimensions 
10. Friction 


A one-half hour exam was given on each of these ten 
topics, with the student being required to pass all ten 
of these basic exams in order to obtain a C in the 
course. Failing to pass one of the basic exams resulted 
in a D grade; failing to pass two resulted in not passing 
the course. In addition, supplementary exams were 
given on wrenches and belt friction, and advanced exams 
were given in the areas of centroids and centers of 
mass, frames and machines, trusses, beams, equilibrium 
in three dimensions, and friction. After satisfying 

the requirements for a C, a student could obtain either 
a B or an A by passing two or four supplementary or 
advanced exams, respectively. During the quarter the 
student was given twenty opportunities to take exams. 
The classes met once a week, during which time some of 
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the more important or more difficult concepts were 
discussed and techniques to solve various types of 
problems were demonstrated, Tests that had been failed 
could be repeated. Finally, all students were required 
to proceed in the course at a minimum pace, which was 
established by the weekly lectures. 


Many of the ills present in the 1971 system were 
corrected by the modifications made during the 1972 
school year. An increasing number of students finished 
the course early, and the grades in the programmed 
courses became substantially higher than in the regular 
sections. However, the complaint is frequently heard 
that the programmed sections required more outside 
preparation from the student. With a little reflection 
this is to be expected, since less time is spent in the 
classroom, 


During the winter quarter of 1973 a modification was 
made in that the exams were no longer graded on a pass~ 
fail basis, and advanced and supplementary exams were 
discontinued. There had been an observed tendency by 
the student to do enough on the exams to just pass 

them, since there was really no reward for doing any 
better. Now, after the student has satisfied the 
requirements for a C, his grades are averaged and grades 
of A and B are determined by his overall average. 


Tentative plans have been made to program nearly all of 
the statics courses at Georgia Tech in the near future. 
Somewhere between 800 and 1,000 students will probably 
take statics by this method during the calendar year. 
However, usage of the two-volume text during the current 
eee has shown that some minor revisions are still in 
order. 


In conclusion, programmed instruction at Georgia Tech 
in statics and dynamics has brought to the student a 
unique and new experience. He now has at his immediate 
disposal a course that he can readily tailor to his own 
desired rate, and use to challenge his ability and 
interest at any level that he desires. Further, the 
instructor is freed of many of his routine teaching 
duties and can now use his time more advantageously 
helping those students that need help on an individual 
basis. Thus, progremmed instruction at Georgia Tech 
has helped students to more readily relate their study 
to the achievement of the course objectives. 
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APPLICATIONS OF MATRIX FUNCTIONS 


J. S. Frame 
Michigan State University 


FUNCTIONS OF A DIAGONABLE MATRIX 

A large class of problems in applied mathematics and 
engineering can be modeled by linear algebraic, 
differential or difference equations whose solutions 
can be expressed compactly and efficiently by the use 
of analytic functions of a matrix. Computer solutions 
can then be programmed using certain algorithms and 
techniques associated with matrix functions. 

For example, the linear differential system 

dX/dt = AX, X(0) =C (Eien) 
A, a constant n x n matrix, has the formal solution 
X= eAtc, and the system 

d2x/dt? + B2x = 0, X(0) = C, (dX/dt), = V (1.2) 
has the formal solution 

X = (cos Bt)¢ + (BU sin Bt)v. (1.3) 


Here the functions of matrices A and B* may be 
defined by series 


eto r+ at+ Ate2s2z1+ ++ - (4) 
cos Bt = I - B*t2/21 + B4e4yg) - + + - (Td 
Bol sin Bt = It - B2t3/31 + B4tSysi- - - = (1.6 ) 


so that only positive integral powers of A and B2 

are involved in the solution. There are important 
computational advantages, however, in expressing these 
and other analytic matrix functions (defined by poly- 
nomials or convergent series in A or A~+), in finite 
form as polynomials of degree less than n in the given 
constant matrices. This is easily done if the given 
matrix is a diagonal matrix A, with diagonal entries 
4 called its eigenvalues. For the equations 
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for any positive integers k, and also for negative 
integers if |A| # 0, enable us to sum the infinite 
series for f(\;), when they coverge at each \j and 
write 3 ; 


Alt 
emt) "6%. 88) £41) 0 +250 
Ajt 
2 
i°} e Dias) 10) E(Ag)*+* 0 
elt = f(A) = a 
xt 
0 Oss eh 0 O ++-£(An)}. 


Suppose that A has s distinct eigenvalues j, 42, 
+++, As, where A; appears with multiplicity mj and 
=mj =n. Then the Lagrange interpolating polynomial 
defined by 


aj) = -a)/@5-ay) 


hj (1.9 ) 
has the property that q;(A,) = 6;,. Hence, the matrix 
function : 4 
aj (4) 0 
By = aj(A) = aj Aj) (1.10) 
0 95 On) 


is a diagonal matrix having m; 1's in the entries 
occupied by .’; in A, and org elsewhere. We find 


j 
s 
= Ej =I and (ale boy 
jel 
By+B, = Ej 85x, (1.12) 


so the E; are called mutually orthogonal idempotent 
matrices, Since the Ej are polynomials in A that 
satisfy the equations 
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show that for A= Aj = 4, the equation (4I-A)S; = 0 
has one solution vector, whereas for A = A2 = 1, the 
rank of I-A is 1, and there are 3-1 =2 independent 
solutions Sz, S such that [I-A][S2, S3] = 0. The 
modal equation Rs = SA becomes 


Ze SEP SE Re zo 1 z 0} |}4 0 0 
Wore a *2' +1 2 Se Dstt LPO Leroy Tote sa. 
aD. aa tema) oe SLO Oped 
(1.29) 


A normal matrix is an nxn (real or complex) matrix A 
such that 


A*tA = AA*, (1.30) 


where A* denotes the complex conjugate of the transpose 
of A. The class of normal matrices includes hermitian 
(including real symmetric) matrices, such that A* = A, 
skew hermitian (including real skew) matrices, such that 
A* = -A, and unitary (including real orthogonal) matrices, 
such that A* = A-l. Their eigenvalues satisfy respec- 
tively 


Xj = Aj real, for hermitian matrices, (1.31) 


Xj = “Aj pure imaginary, for skew hermitian matrices, 
and Xj = aj of absolute value 1, for unitary matrices. 


Every normal matrix is diagonable by a unitary modal 
matrix S, from whose columns we can form the mutually 
orthogonal idempotents $jSj* of rank 1. For a root 

aj of D(A) = O of multiplicity mj> 1, we must add 

mj such idempotents to obtain the Constituent idempotent 
Aj. If Sj; is an eigenvector of A, but is not neces- 
sarily of dength 1, the idempotent is 


Ag = ES, S,*/(S_*S,), summed for Ay = Aje (1532) 
Given a diagonal matrix Aj we can_expand certain known 
functions £(A) = 1, A, 22,°++-, aS“l in terms of the 


unknown idempotents, and then solve for Aj, Az,<¢*s 5 As 
in terms of the known powers I, A, AZ, +++, AS"1. In 
the example above with 1 = 4, 2 = 1, and A in (1.26) 
we have 


for £(A)=1:1=+1+A,+1+ A2, and (1.33) 
for £(A)=A’A:A=4Az,+1- A2. (1.34) 
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Hence, 
tre ek 2-1 -1 
A-I 1 A-41 .1 
Ay = AT = 3d a Ajy Ages TF 3 Fl 2 et) (1.35) 
BT ee | AS = 2 


APPLICATIONS TO LINEAR STATE MODELS 


A common problem which arises in network theory and in 
more general types of linear systems is the solution of 
the equation 


dX _ (2.1) 
Set AX + BU(t), 


where X is an unknown n x 1 vector function of time, 
U(t) is a given mx 1 "control" vector function of 
time, and A and B are constant n-rowed matrices of 
n and m columns, respectively. 


If A is diagonable with modal matrix S and diagonal 
spectral matrix A, then the change of variables X = SY 
reduces equation (2.1) to 

sat = asy + BU(t). (2.2) 


Left multiplication by R= S"! yields the simplified 
system 


SY _ ay + RBU(t) or i = GY, + RzBU(t) (2.3) 
dt at i j . : 


Thus, the variables are separated, and the explicit 
solution is 


Kathy Rett} 
y,(t) =e) y,(0) + fre) “Rj BU(t)dr. (2.4) 
j j ° J 


As t increases, the leading term, obtained for B= 0, 
will tend to 0, oscillate, or become infinite, accordingly, 
as the real part of the eigenvalue j is negative, 

zero, or positive. The linear system (2.1) is called 
stable if none of the eigenvalues ,; of A have positive 
real part, and no multiple eigenvalte is purely imaginary. 


The solution to (2.1) can be expressed directly in terms 
of A, using functions of A, in the form 


18 


x = eAtx(o) + fteA(t-")pucryar, (2.5) 
ce) 


if formula (1.20) is used to express the function eAt 


Ast 
in terms of the function values e J and the idem- 


potents Aj = SjRj- 


APPLICATIONS TO MECHANICAL VIBRATIONS 


The differential system 


2 
d“X dx 
oe a xO) =¢, 71,5 ¥ (3.1) 


furnishes a mathematical model for a variety of 
mechanical vibration problems, such as the vibration of 
airplane wings or of multistoried buildings, where the 
coordinates xj may represent displacements of a 


sequence of points that are equally spaced at equilibrium. 


A simple example concerns the longitudinal oscillations 
of n equal masses attached to a uniform horizontally 
stretched spring of length L = (ntl)h, at distances 
Kh from one of the fixed ends when the spring is at 
rest. If xj denotes the subsequent horizontal dis- 
placement of the ith mass from equilibrium, and k‘¢m 
is the spring constant, then by Hooke's Law we have 


=) 3 Or 0 


2 
md Ear Nid AWE 2 te didisdite Oil “eK (3.2) 
dt? 


This equation may be rewritten in the form 


ax 


peeks k?[21-A]X = 0 (3.3) 


and be solved by (1.3), where B = k2(2I-A), and A is 
the matrix (3.5) with entries 1 immediately above and 
below the diagonal and 0's elsewhere. Now the 
trigonometric identity 
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sin (i-1) j® + sin (i=1)je = (2 cos j@) sin i(j@) (3.4) 


shows that A Sj and Sjrj will have the same ith 
entry if 


0 10 0+-+0 sin je 
yy One Sh 10) 095.0 sin2je 
A= a 
010 10 J | sin3je 
ar GB veh Masia 
‘0, p> oie 10 sin n je 
; © 
Aj = 2 cos jo, @= Tye (3.5) 


Since A is real symmetric, it is normal, and its con- 
stituent idempotents Aj have the form SjSj*/(Sj*S5)- 


The squared length of Sj is 


n 2 1 n 
Sj*Sj = © sin*ij@=> f ¢@ -. cos (2ij8)) (3.6) 
i=l 2 i=0 
ney pao bias (nti) ccttive o 
ergs _= cos 2ijn/(n+1) =o . 
i=0 


Thus, A has the constituent idempotents 


2 
Ajsaey (8985"° tSi7) 


To compute the eigenvalues vj of B= k(2t-a)}/2,, we 
replace A by its eigenvalue” Aj = 2 cos j®, and find 
uj = k(2-2 cos je)1/2 = 2k sin jg, (3.8) 
where @= 0/2 = 1/2(n+1). 


Finally, using (1.3), we obtain the explicit solution 
OF) (Bed) 2 


2k sinjg n+l 


n 7 Soir 28; 
X = 2 [(Sj*C)cos(2kt sin jg) + (s;*v)Sin(2ktsinj eas 
an (3.9) 
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If the more general equation (3.1) is to be studied 

for large values of n, it may suffice for practical 
purposes to find only the lowest oscillation frequencies 
Wl, W2, ‘** With their corresponding modal vectors Sj- 


A vector ets; satisfies the differential equations 
(3.1) if and only if 


(-w2 M+ K) Sj = 0, or (we? 1 - K-12 M) Sj = 0. (3.10) 


The numbers dy = we? represent for k = 1,2,3 the 
largest eigenvalues of the matrix A = K~4M, for which 
Sk is a corresp ing eigenvector. Since e symmetric 
matrix mien Pag the eigenvector maPas,® it 

has the corresponding idempotent 


1/25. g.*yl. sMS; 
wi/2s; sj*m1/2/(s;Ms;). 
Thus, the idempotents of A are 


$585 *M/(S3*MS;). (3.11) 


If the largest eigenvalue dy has multiplicity 1, 

we can calculate it without finding and solving the 
characteristic equation by using the following 
iteration procedure: Starting with the first column I; 
of the unit matrix, we successively multiply by A 

and divide by the first non-zero coefficient. Then 
the limiting vector is the eigenvector S,, and the 
limiting divisor is the largest eigenvalie a: 

From (1.20) we have 


(A/aq)¥1y = [Ay + (A2/A1)KAg + +++ + Cin/dq)KAQT2, (3.12) 


and for large k the right hand side approaches Ailj, 
which is an eigenvector that is unaffected by 
multiplication by A/Aj. 


A DIFFUSION PROBLEM IN FLUID FLOW 


In a certain chemical engineering process a fluid flows 
continuously with constant velocity C in the positive 
s_ direction through a long cylindrical tube or reactor 
of length L with constant circular cross-section. At 
time t= 0, a new fluid is introduced at the input end 
and flows through with the same velocity. Part of the 
old fluid is washed out, but part diffuses back through 
the new fluid. The problem is to determine the concen- 
tration x(s,t) of the new fluid as a function of 
positive s and time t. 
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We first assume that the variation of x can be 
described by the Taylor diffusion model 


2 
aX = -¢ 3X+ pox 4.1 
at us 982 G-9) 


with velocity C and dispersion coefficient D, subject 
to the initial and boundary conditions 


x(s,0) = 1 for O<s<L, and (4.2) 
x(0,t) = 0, ax(L,t) = 0, for t>0. (4.3) 
2 
at 


Next we approximate the partial derivatives with respect 
to s by finite differences, taking s = kh, h = L/n, 
and obtain a system of first-order ordinary differential 
equations 


dxy Xk+1 ~ Xk-1 X41 ~ 2x * Xk-1) 
wc se +D = (4.4) 


where x,(t) = x(kh,t) = 1 for t = 0, and 


Xo = ‘OF Xnt1 72x) +°X,-7 20 for t>0. (4.5) 


With a new time variable + and parameter p 
defined by 
t = tD/L2, p = Ch/2D = sin q, (4.6) 


the differential equation (4.4) becomes 


dx 
n 7k = (1#p) x21 - 2x_ + C-P) yA. + (4.7) 
We eliminate x and x,,, by the boundary conditions 
(4.5), introduce the n-vector X with components 
Xy, °**, X) and combine equations (4.5) into a 
single matrix equation: 


72 dX LL 
no’ $t = - MX. (4.8) 
The tridiagonal coefficient matrix -M is not symmetric, 


but may be made symmetric by transforming by a diagonal 
matrix B. We define 
a= (i1-p2)!/2 = cos g, b = (1-p)/q = tan(n/4-q/2), 


© ean 4yi/ 2, (4.9) 
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M = 21 - q BUAB, B = diagyb, b2+-+, b®-1, wA/cy, (4.10) 


where 
cay aro) 0 0 x] 1 
“Soot jon ila ha aba) X2 1 
Aj }O0 1 O70 0 X= x3 Xs 1 (4.11) 
Owe ka Ott ic 2 fe 
ONOve 2s Ov c. 2b x 1}. 


The required solution of the differential equation that 
satisfies the initial and boundary conditions may be 
written 


2 
coe Ge a (4.12) 


To the eigenvalue A. = 2 cos 6; of the symmetric 
matrix, A corresponds the eigenvalue 


uj = 2-qhj = 2(1-cos @ cos 65) (4.13) 


of M. The eigenvector Sj; of A in (3.5) must be 
modified slightly to give’ the following eigenvector 
for A in (4.11): 


* = [si +, si ESNet in(n-1)e- i te 

8; [sin 6j, sin 285, » sin(n 1)0;, (sin n Aa 
where 6; is not jn/(nt1) as before, but is 
determinéd by an equation derived from the last row of 
A, namely 


© sin(n-1)8j + 2(b/c) sin nos = (2 cos 6; sin n 93)/c 
j j j J 


(4.15) 
which may be rewritten as 

sin(n+1)e; - 2b sin n 0; + b? sin (n-1)0j = 0. (4.16) 
The constituent matrices for A and M are 

Aj = SjSj*/(S;*8j) for A, Bo1A5B for M, (4.17) 


where by the aid of some complicated trigonometric 
identities we find 
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31), and (4.18) 


85*Sj = (n/2)(1 + 5 


vj = 03 n/2p = nj n? D/CL. (4.19) 


Restoring the time variable t from (4.6) in the 
solution (4.12), and using vj; from (4.19), we obtain 
on expanding the exponential ~ function 


n ~vjCt/L 
X= 2 rj e Bls;, (4.20) 


where vy denotes the scalar quantity 


Tj = Sj* B Xo/(Sj*S;), (4.21) 


which can be reduced by a series of trigonometric 
identities to 


rj = (sin 95)/ (tan ) (1 + vj)- (4.22) 


To obtain A; = 2 cos 6;, and hence yj and v3, we must 
solve the characteristi¢ equation (4.16), we observe that 
the left side 9§ farie), is the imaginary part of 

2n(zl1/2 - bz-1/2)2 for z = e18. Noting that 

cot (y/2) = (1+b)/(1-b), we can express condition (4.16) 
in the form 1 

jx =n 6j + 2 tan™”(cot(@/2)tan(0j/2)). (4.23) 


We determine q@ = sin-1(ch/2D) from (4.6), solve for 

6j by iteration in (4.23), calculate yj = (1-cosacosaj)- 
n/sinw, find rj from (4.22) and compute the solution 

X from (4.20). Then 


n 7vjCt/L k 
X= 2 e b°% sin k 0j sin 6j/(l=vj)tan cp, 1<k<n. 
bie % ACh * (F284) 
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RESEARCH SYNDROMES 


S. B. Batdorf 
The Aerospace Corporation 
El Segundo, California 


I would like to say a few words about the quality of 
mechanics research, or more precisely of engineering 
research in general, starting with a true story of a 
technical meeting I attended about S years ago, This 
meeting was set up for the purpose of acquainting the 
scientific community with the latest results in upper 
atmosphere bomb phenomenology. It was attended by about 
a thousand people and lasted 3 days. Each day there 
were about 15 papers, and each paper employed 15 to 20 
figures or tables, That is nearly a thousand figures 
and tables - almost entirely raw data, either experi- 
mental results or computer calculations. Nobody was 
allowed to take any notes for security reasons. When 
the meeting was over I don't believe any of us could 
have talked more than 10 minutes on what we had learned 
in the 3 days, which is a damn poor rate of return on 
investment. 


On the morning of the second day, a paper was given on 
a computer analysis of a 1 megaton burst occurring at 
30,000 foot altitude, There were scads of graphs giving 
the results, and I found myself struggling to make . 
sense out of them. One, for instance, was a logarithmic 
plot of the time after burst at which the shock wave 
would arrive at various distances from the point of 
burst. I found myself using one finger to help estimate 
the slope of the line to see if it was 2/5, as predicted 
by Taylor's simple theory of point explosions published 
in the Proceedings of the Royal Society about 20 years 
ago. But before IT could quite manage it, the slide 
disappeared, and we were looking at the next one. 


It turned out that in the break following this paper, I 
wound up in the coffee line behind the sponsoring . 
agency's chief scientist. I mentioned that we were being 
inundated by a flood of undigested data that nobody could 
assimilate and take away with him. I said that as soon 
as I got my copy of the proceedings I would get out a 
ruler, check slopes, sort out which results agreed with 
what I expected, which didn't, and try to interpret the 
deviations. But instead of one thousand individuals each 
doing this for himself, I said, why doesn't the author 
do it? He said he hadn't thought about it before, but 

it seemed like a good idea. He would try to see that 
this is done in future meetings. 
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The very next paper was one dealing with the problem of 
an astronaut in a space capsule at 100 miles altitude 
attacked by a 1 megaton co-altitude nuclear bomb. The 
author calculated the neutron flux through the astronaut's 
body for various burst distances. He told us with some 
pride how he used FORTRAN to determine the neutron flux 
inside the space capsule, and he gave the results as a 
logarithmic plot of flux versus distance with about 50 
computed points. They fell on a very good straight line, 
I checked the slope and found that it was minus 2. I 

then noted that the chief scientist was in the row 

behind me and moved over to speak to him. "Look," I said, 
"I believe our speaker has just proved the inverse square 
law!"" He look surprised and then burst into a laugh 
which must have startled the rest of the audience. I 
suspect that less than 1% of those there arrived at this 
or any other interpretation. The reason is very simple. 
We are so accustomed to barrages of meaningless data 

that nearly all of us turn off our minds and stop thinking. 


I chose this example because I think it illustrates rather 
well a number of syndromes which I believe afflict 
engineering research today, and I would like to spend 


the rest of my time identifying and discussing some of 
these. 


One syndrome is the standard practice of writing reports 
designed to impress rather than inform. In the case of 
the astronaut and the space vehicle, the most meaningful 
way to report the results, and one that we could have 
carried away with us, would be to say that for the 
particular space capsule the fraction of neutrons that 
got through and hit the astronaut was x$. When you 
come right down to it, the shielding factor (which 
actually wasn't mentioned) was really the entire useful 
content of what the investigator did. But that doesn't 
sound nearly as impressive as to report a computer 
program that probably employed a Monte Carlo approach to 
calculate the shielding, and demonstrate that it had 
been successfully employed in 40 or 50 cases. 


A second syndrome which is related to the first is our 
habit of suppressing our intuitive processes, or at 

least destroying the evidence. A former president of 

the American Physical Society, named W. G. F. Swann, 

once wrote a book called the Architecture of the Universe, 
and I would like to read an excerpt from that. A 
scientific researcher can be likened to two distinct 
animals, Ahead there travels a dog with keen scent and 
not over critical capacity who pokes his nose into every 
bit of intellectual garbage he can find on the wayside. 
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He ferrets out from every conceivable place new morsels 
which no dog has ever chewed before and gnaws out as 
much of their content as he can, And behind the dog 
there walks, clad in the most conventional apparel, a 
calm and dignified gentlemen who views very critically 
everything the dog picks up, and makes him put it down 
again if it is found wanting, but pockets it in his 
immaculate coat if it is of worth. It is not to be 
forgotten however, that is is usually the dog that finds 
the bone. The dog's name is Intuition and his master's 
name is Precise Thought."' In writing papers we generally 
eliminate all traces of the intuitive processes which 
lead to an important discovery, and give only formalized 
results. Intuition is somehow considered not quite 
respectable, possibly because when properly employed it 
often makes a laborious formal demonstration of something 
seem totally unnecessary. 


Incidentally, in his speech following his acceptance of 
the Timoshenko medal a few years ago, Bill Prager 
reminisced about the old days in Gottingen and the con- 
trasts between two of the leading figures there, Prandtl 
and von Mises. Prandtl always looked at things phse realy 
and sought to make things intuitively evident, while 

von Mises always insisted on rigorous mathematical 
analysis before he would accept any results. Both 
approaches are of course important to progress, but Prager 
expressed the view that today the pendulum has swung too 
far in the direction of formalism. Inventiveness is 
greatly enhanced by a careful cultivation of a broad 
intuitive grasp, corrected whenever it is in error. And 
for whatever it is worth, I suspect that in the annuals 
of fluid mechanics Prandt1's name will be more illustrious 
then that of von Mises. 


Another syndrome is our tendency to divorce our technical 
thinking from our everyday experience, As a result, we 
lose perceptiveness and suffer some degradation in the 
valuable faculty known as common sense, at least as 
applied to technology. As a simple example, take Kelvin's 
theorem about circulation. It is sometimes expressed as 
follows: Vorticity cannot be induced in a perfect fluid 
by forces derivable from a potential. Yet everybody who 
has ever poured tea in a cup or flushed a toilet knows 
that a whirlpool is set up. And everybody who knows 
Stokes! theorem relating the surface integral of the 

curl of a vector to the line integral around the surface 
knows that a whirlpool implies the existence of vorticity 
somewhere inside, But in most cases experts in fluid 
mechanics, to whom I have posed this paradox, had never 
previously been aware of any contradiction between their 
accepted technology and this everyday experience. And 
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in explanation, 90% have adopted the view that in teacups 
and toilets circulation does indeed exist, but it is 
either because water isn't a perfect fluid or because the 
forces applied aren't derivable from a potential, I 
think both explanations are wrong, but I'll leave it to 
you to come up with the correct one. 


Then there are the cases where people apparently don't 
really want to clean up a field--that would deprive them 
a pleasant and secure job. Many of you have probably 
heard the story about a young lady who went into a depart- 
ment store to buy some material for a nightgown. She 
wanted something very soft and gauzy and the saleslady 
finally found just what she wanted. When asked how much 
was needed, she requested 20 yards. The saleslady was 
aghast and said, "My dear young lady, you obviously don't 
know much about making clothes, All you need for a 
nightgown is 2 1/2 yards." The sweet young thing replied, 
"But you don't know my husband. He's a research and 
development engineer, and he would much rather look for 
something than find it!" 


A very serious syndrome is our habit of turning our 
thinking over to the computer. I am personally of the 
opinion that if the computer had reached its present 
stage of development in the time of Sir Isaac Newton, 
the science of mechanics as we know it today would not 
exist. The computer has such an enormous capacity for 
incorporating unimportant details as well as the main 
features of a system into an analysis, that the art of 
idealizing things into particles, rigid bodies, weightless 
strings, frictionless pulleys, and all the rest might 
never have developed; nor might the concepts of kinetic 
energy, potential energy, momentum and all the other 
abstractions which shape our thinking in mechanics. 


An amusing example of this habit of relying entirely on 
computers occurred in the early days of the space program. 
When the second moon probe was launched (the first blew 
up), there was a great deal of excitement and newscasters 
gave reports every quarter or half hour on the probe's 
altitude and velocity. For several hours it was taken for 
granted that the probe would reach the moon, the only 
question being whether it would hit it or miss it. Being 
somewhat skeptically inclined, however, I took out my 
slide rule when I heard the first report on altitude and 
velocity about fifteeen minutes after launch, calculated 
the potential and kinetic energies of the probe, and 
calculated the maximum altitude that it could reach, I 
found that it couldn't go even half way. For hours after 
that, newscasters kept reporting the probe's position 

and velocity with no doubts expressed about its ability 
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to reach the goal. But eventually the truth was out, 

The probe reached 80,000 miles, 1/3 of the distance to 
the moon, returned, and fell in the ocean. This type of 
calculation takes less than 10 minutes, and it was 

amazing to me that when so many people in the United 
States were perfectly capable of doing it, nobody bothered; 
or if some did, they didn't communicate their results to 
the news media. 


Another common syndrome is losing track of the fact that 
there is no such thing as an exact solution to a real 
physical or engineering problem. The only thing that 
permits an exact solution is some idealization of the 
physical or engineering problem, and from a practical 
point of view there is no great advantage in reducing the 
error in the solution of the idealization very much below 
the error by which the idealization fails to represent 
the problem itself accurately. There is really nothing 
new in this observation. In fact, I believe it was 
Aristotle who once said that it is the mark of an 
educated man not to expect greater accuracy than the 
subject permits. We characteristically lose sight of this 
advice, however, and carry out the idealization to as 
many significant figures as we can manage--very often 

at excessive cost and with unnecessary time delays. 


After some experience with a given idealization, the 
idealization may assume greater importance in our minds 
than the problem itself. This reminds me of the story 

of the policeman who was on his beat early one morning 
about 1:15, and he saw a man in a tuxedo crawling around 
on his hands and knees under a street light. He went over 
thinking the man must be drunk and that he ought to run 
him in, but discovered it was a man he had known for 

ears and rather liked. Upon inquiry, it turned out that 
he had lost his house keys. Being of a kindly disposition, 
the policeman got down on his hands and knees to help 

look for them. After about 10 minutes he became dis- 
couraged and asked the man where he lost them, "Oh, 
the man replied, "over the across the street." The 
policemen asked, "Then why don't you look over there?" 

to which the inebriate replied, "Don't be silly. It's 
dark over there!" I sometimes think that's the story of 
research. We tend to work on the things we feel we under- 
stand instead of the problem that really needs a solution. 


Part of our excessive dependence on computers may come 
from an instinctive desire for good accuracy. But the 
accuracy that is needed depends on the problem, and 
often in the initial stages of an investigation very 
crude numbers are quite adequate. For instance, in the 
early days of the communications satellite program, some 
inventive people at Lincoln Laboratory and Rand 
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Corporation suggested that a very effective, practically 
indestructible passive communications satellite system 
allowing instant communication with other continents 
would be to have a cloud of tiny dipoles traveling in 
orbit and filling a complete circular path around the 
earth something like Saturn's rings. In order to get 
as many dipole reflectors into orbit as possible subject 
to a given weight limitation, the individual elements 
must be very thin. As a result, the radiation pressure 
of the sun becomes important in perturbing the orbital 
path. 


At that time I was serving in the Pentagon as Technical 
Director of the Military Communication Satellites Program, 
and was very anxious to get an answer to the longevity 
of the system. On a Thursday afternoon, I remember 
calling some people who were calculating the orbital 
perturbation, and learned that they were putting the 
finishing touches on their computer program and expecting 
some answers the following Tuesday, Unwilling to wait, 
I idealized the problem as follows: The largest pertur- 
bation will come when the sun is in the orbital plane 

of the needles. When that happens we can divide the 
nominally circular orbit into four arcs. In one of them 
a needle is in the earth's shadow and is thus unperturbed. 
In the corresponding arc on the sunny side the average 
velocity of the needle is at right angles to the line 
joining the sun and the earth. The two remaining arcs 
consist of one in which the average direction is directly 
away from the sun. We can thus approximate the perturbing 
effect of the earth's sun radiation pressure by calcu- 
lating the total impulse that will be given as a single 
pulse at the center of each of the arcs. It is easy to 
show that each of the three impulses has the effect of 
changing the orbit to an ellipse with its major axis at 
right angles to the line of joining the earth and sun, 
and that the ellipticity increases with each pulse. Thus 
the apogee increases and the perigee decreases until 
finally the orbit dips down into the earth’s atmosphere, 
after which, of course, the remaining lifetime is very 
short. I worked out a formula for this on Friday after- 
noon, and was pleased to find the following Tuesday that 
the computer verified my results. Here was a case in 
which the inaccuracy involved in replacing continuous 
pressure by three impulses was far less than the uncer- 
tainty in the sun's radiation pressure on each dipole 

due to its random orientation. 


By now many in the audience will probably feel I am 

anticomputer, Actually I'm not. The computer is of 
inestimable value in solving many important types of 
problems. For some it is even indispensible, But it 
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should be kept in its proper piace and must not be 
allowed to replace physical thinking, intuitive insight 
and a broad grasp of fundamentals. There is a consider- 
able threat to these arising from analytical solutions 
also. When solutions are presented as rans intergrals 
or complicated functions of variables which have been 

so transformed that you can't see clearly how they are 
related to the actual physical parameters, the results 
are little better than cooking recipes. Like computers, 
such solutions have a place, but they leave the essential 
job of conveying understanding still to be accomplished. 


What I do think is that a complicated and expensive solu- 
tion should only be resorted to when the situation 
warrants it. Further I believe that each such solution 
should be preceded whenever possible by a simple solution 
which scopes the problem, conveys a basic understanding 
of its main features, and shows whether and for what 
values of the parameters a more accurate solution is 
worthwhile. There are some people who have developed 

a remarkable capacity for stripping a problem down to 

its essential elements and getting simple approximate 
solutions. Most of the people I know on the President's 
Scientific Advisory Committee and high level defense 
technical councils are back-of-the envelope artists in 
their own right. I don't believe this is a mere 
coincidence; rather, there is a cause-and-effect relation- 
ship between technical competence and breadth on the one 
hand, and the practice of viewing a problem in its 
simplest terms on the other. 


This is far from a complete discussion of research 
syndromes, but I have gone on long enough. If we could 
lick those I've mentioned, I believe it would improve our 
performance remarkably. So let me just close by 
summarizing. 


We should carefully cultivate our basic physical under- 
standing of the phenomena we deal with, and recognize the 
important role of intuition in discovery. Technical 
papers should be written so as to inform rather than 
impress. We should integrate our daily experience and 
our technical thinking for the additional insight and 
practical common sense that can result. The technical 
approach employed should be properly tailored to the 
needs of the occasion, using simple approximate methods 
to scope a problem and achieve a basic understanding, 
followed by relatively complex and expensive solutions 
whenever the situation warrants it. Finally, we should 
be guided in our choice of problems by their importance, 
whether theoretical or practical, and not be shackled by 
either our past history or tradition. If we can really 
follow all of these admonitions, I believe it might 
possibly result in something of a renaissance in mechanics. 
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APPLICATION OF MATLAN TO STRUCTURAL ANALYSIS 


M. A. Eisenberg and G. W. Hemp 
University of Florida, Gainesville 


INTRODUCTION 


It has long been recognized that the computer is an 
indispensible tool in the stress and deformation analysis 
of all but the most simple structural configurations. 
Continued development of digital computers has made the 
current generation faster, more reliable, and more 
economical to operate than ever before, However, the 
primary difficulty still remains--communication with 

the computer. In recent years increased attention has 
been given to the development of problem-oriented and 
procedural languages in order to alleviate this difficulty 
(Fenves 1969; Murphree, Melin and Fenves 1966; Melin 
1966). In particular, problem-oriented languages (Hart 
and Hargraves 1960; Miller 1961; Fenves et al. 1964) 
simplify greatly the task of the programmer. However, 
this simplification is achieved by specialization and 
limits the number of workers to whom a particular 
language would be of use; this requires a given worker 
to be familiar with many computer languages in order to 
handle the full scope of EF eeRlane with which he is 
concerned. 


Thus, high level procedural languages of a general 
nature such as FORTRAN IV and PL/1 must continue to be 
the primary tool of most computer users. However, these 
languages should be supplemented by others which are 
designed to handle a selected class of operations 
conveniently. One such language is MATLAN (Matrix 
Language), which was developed by IBM Germany for use 
with 0S/360 and was released in December 1968.* The 
features of this language are such that it should prove 
to be extremely useful as a classroom and research tool 
in the study of matrix methods of structural analysis. 
The facility with which structural mechanics codes may 
be generated by students or researchers with limited 
programming experience will be demonstrated below. 


MATRIX LANGUAGE 


Basically, MATLAN is a programming system which is 
designed to simplify handling of and computation with 
matrices. The language contains provision for matrix 


¥System/300 Matrix Language (MATLAN) (360A-CM-05X), 
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generation, manipulation, arithmetic and program control, 
and may call MATLAN or FORTRAN subprograms. There are 
two primary features which make the language particularly 
attractive. The first is the extreme ease and conciseness 
with which matrix operations may be coded; this feature 
will be demonstrated in an example to follow, The 
second major advantage is that the automatic storage 
management scheme controls the dynamic flow of data within 
the system. It is responsible for dynamic allocation, 
renaming and canceling of matrices both in core and in a 
random access file. Matrices are automatically stored in 
array or coordinate form according to the density of data 
and are dynamically allocated to core or direct access 
file during execution. Two versions exist for all 
algorithms. One version, the core algorithm, occurs when 
sufficient core storage is available for the arrays; 
otherwise a segmenting algorithm is executed. Thus the 
matrix operations are not only simple to code, but are 
executed in a more efficient way than could be expected 
in a FORTRAN code written by one who, though competent 
in matrix concepts and in the application area, is a 
novice or occasional programmer, 


STRUCTURAL ANALYSIS PROBLEM 


In order to illustrate the application of MATLAN to 
structural analysis problems, a general method of elastic 
analysis of redundant structures by a matrix force method 
due to Wehle and Lansing (1952) has been coded. The 
resulting program is used to analyze an idealized wing 
structure. 


Let the structure under consideration be divided into 
easily analyzed members, and let the loads acting on the 
individual members be given by {Q} = {Q'} + {Q"}, where 
{Q'} represents s statically determinate loads and 

{Q"} represents r redundant loads. Then s _ independent 
linear algebraic equations of equilibrium may be written 
for the structure, These equations may be solved for 

{Q'} and combined with r trivial relations for the 
redundant loads to give the s + r equations 


nie 
=e ty, a) 
where the P are known applied loads. 


The redundant loads are determined by minimizing the 
strain energy 


U = FIQUIAN(Q) = FIP} Quy(B}(PL Qn]? (2) 
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where [A] is the influence coefficient matrix and 


‘ 
Bit Be 1207 ae 
[B] = Big > [616 ] {AJ[GiG]. (3) 
The energy is extremized by requiring the vanishing of 
aU/aQ". The resulting equation in terms of {P} and{Q"} 
may be solved for {Q"} to obtain 
4 
{Q"} =-[B']7*[B°]{P} = [c]{P). (4) 
From (1) and (4) 
{Q} = [G’ + G” C]{P}. (5) 
To complete the solution, Castigliano's theorem is 
applied to obtain the displacements of each of the n 
applied loads P. Upon differentiating U with respect 


to the applied loads and making use of the solution 
for Q" one finds 


{D} = [BB] [Hm. (6) 
Thus, given [A] and [G], the internal load distribution 
and deflections of the applied loads may be computed for 
any structure. 


MATLAN CODE 


Hee following code will implement the analysis described 
above: 


1, MAIN WEHLE 

2. READ (A,G) 

3. HEADING "STATIC INDETERMINATE STRUCTURE MATRIX 
SOLUTION! 

4, TEXT "SUM OF MEMBER FLEXIBILITIES', C 

5. WRITE A 

6. NEWPAGE 

7. TEXT ‘UNIT LOAD DISTRIBUTION’, C 

8. WRITE G 

9. NEWPAGE 

10. MULT A,G,X 

11. TRANS G,GT 

12. MULT GT,X,B 

13. EXSUBM B, (1,1), (4,4),B1 

14. EXSUBM B,(1,5),(4,3),B2 

15. EXSUBM B, (5,1), (3,4) ,B3 

16. EXSUBM B, (5,5), (3,3), B4 

Lat DIV: B4,B3,Y 

18. COPYC Nye 

19. EXSUBM G, (1,1), (17,4) ,G1 
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20. EXSUBM G, (1,5), (17,3) ,G2 

21. MULT G2,C,2 

22, ADD G1,Z,LOADS 

23, MULT B2,C,Z 

24. ADD B1,Z, INFCOF 

25, TEXT "REDUNDANT STRUCTURE UNIT LOAD 
DISTRIBUTION" 

26. WRITE LOADS 

27, NEWPAGE 

28, TEXT "INFLUENCE COEFFICIENTS! 

29. WRITE INFCOF 

30. END 


In this program, the actual computations are described 
in statements 10-24, the remaining statements being 
devoted to input/output bookkeeping. Statements 10-12 
compute the product indicated in Equation (3). 

Statements 13-16 and 19-20 partition [B] and [G] as 
indicated in Equation (3). The numerical parameters 
which appear in these statements correspond to the matrix 
sizes of the sample problem to be described below. In 

a working program, they would, of course, be treated as 
variables to be read in as data at execution time. 
Statements 17 and 18 perform the matrix inversion, 
multiplication and sign change indicated in Equation (4); 
21 and 22 compute the load distribution according to 
Equation (5); and 23 and 24 compute the influence 
coefficients according to Equation (6). The advantages 
to the student of coding this analysis in MATLAN as 
opposed to FORTRAN or even PL/1 are manifestly evident. 


SAMPLE PROBLEM AND CONCLUSIONS 


In order to test the program, the same example as employed 
in the original Wehle and Lansing (1952) paper was 
selected. A tapered swept-back wing was idealized by a 
conventional sheet-stringer box. This required the 
solution for 17 unknown loads in a thrice redundant 
structure with 4 applied loads. Thus s = 14, r = 3, 

n = 4 and [G] is a 17x17 matrix and [A] is a 17x17 
matrix. The [A] matrix is sparsely populated so the 
MATLAN algorithms will, automatically store it in com- 
pressed mode. Also [B'] can be shown to be symmetric 
and positive definite. These properties could have 
been specified in an ATTRIBUTE statement in the program, 
and the solution algorithm employed in the DIV instruc- 
tion in statement 17 would have automatically employed 
an algorithm appropriate to such a matrix. Since such 
refinements were not incorporated into this sample 
program the solution proceeded as if the equations were 
asymmetric, The computations were carried out on an 
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IBM 360/65 computer running under MFT-II (Multi-program- 
ming with Fixed Tasks). The CPU time expended on the sam- 
ple problem was approximately 7 seconds, including com- 
pile, link edit and execution. 


This example illustrates the ease with which a structural 
analysis program can be written in MATLAN. The student 
never loses sight of the operation which is to be carried 
out by becoming involved in the details of that operation, 
e.g., the inversion of a matrix. With the availability of 
this new language, students without prior programming ex- 
perience learn to apply matrix structural and analytical 
techniques to simple but realistic problems very quickly. 
In addition, the student will have been introduced to a 
useful engineering tool since the language allows one to 
handle extremely large matrices. 
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THE APPLICATION OF ACCELERATION ANALYSIS 
TO THE COMPUTATION OF THE PRIMARY 
FLOW IN TURBOMACHINERY 


Kalman I. Krakow 
Loyola of Montreal 


The absolute acceleration of a particle, ay, in motion 
with respect to a rotating reference frame consists of 
two components: the acceleration of the particle 
relative to its coincident point in the rotating refer- 
ence frame and the absolute acceleration of the coinci- 
dent point. Each of these two components is in turn 
equal to two physically identifiable components. The 
first component is equal to the acceleration of the 
particle relative to the rotating reference frame as if 
the frame were static, a,, and the Coriolis acceleration, 
cor: The second component is equal to the centripetal 
acceleration, 4), and the tangential acceleration, a>. 
Hence, the foll6wing equation is obtained: 


ay = 8% + Acorn * Fcp * Br- 
A C 11] 


This equation is commonly used in mechanics of particles 
and in mechanics of machines. It may also be applied to 
the fluid mechanics of turbomachinery enabling a 
formulation of the governing equations which has the 
following advantages over conventional formulations: 


1. Greater insight into the mechanics of the flow 
is permitted since all terms are physically 
identifiable. 

2. Comprehension of the concepts conveyed by 
equation [1] itself is aided. 

3. Solution of the equations by numerical methods 
may be performed easily. 


Although the writing of a complete program requires too 
much time for incorporation into a fluid mechanics course, 
a student may write one or more of the required sub- 
routines. The algorithm for the solution of the equations 
governing the flow is basically as follows: 


1. A first approximation to the relative streamlines 
is calculated. 

2. Flow angles, streamline curvatures and other 
necessary geometric parameters are calculated, 
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3. The relative total temperature and the relative 
total pressure distribution along streamlines 
are calculated by using the equation of conser- 
vation of energy, the Euler turbine equation, 
and the equation relating temperatures and 
pressures adiabatically. 

4. The continuity equation and the Euler equation 
of motion are solved simultaneously along 
calculation lines, which cross the streamlines, 
to yield a new closer approximation to the 
streamlines and required aerodynamic parameters. 
The format of the Euler equation of motion used 
embodies equation [1] and is the factor which 
makes this analysis unique, 

5. Using the new approximation to the streamlines, 
the foregoing procedure is repeated until the 
repeatability of results is within an acceptable 
tolerance. 


In order to illustrate this application, the following is 
a formulation for two-dimensional radial turbomachines. 
(Similar formulations may be used for axial and mixed 
flow [three-dimensional] turbomachinery, ) 


A rotor is considered to consist of a series of 
channels formed by its blades. The flow in a single 
channel is to be analyzed. (Overall turbomachine 
parameters may be calculated from data obtained from 
the single channel analysis.) The channel is defined 
by the pressure surface of one blade and the suction 
surface of the next blade. In order to facilitate 
relating geometric and aerodynamic data and equations, 
three coordinate systems shown in Fig. 1 are 
required: a Cartesian and a cylindrical coordinate 
system having the center of rotation as the common 
origin, and a natural coordinate system. Calculation 
points, defined by the intersection of the relative 
streamlines (indexed i), and calculation lines of 
constant radius (indexed j), are established within 
the channel. The blade surfaces, i.e., the boundary 
streamlines, are specified numerically by arrays of 
Cartesian coordinates. The coordinates of the con- 
secutive points on each blade surface are taken such 
that they lie on the same calculation line. 


In order to maintain the same equations for both compres- 
sor and turbine geometries, the x component of the velo- 
city is always taken to be in the positive x direction, 
therefore channels representing compressors and turbines 
are to be specified by positive and negative x coordinates 
respectively. Furthermore, the cylindrical coordinate 
systems must be different for compressor and turbine 
geometries, as indicated by the unit vectors i, and a, 
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shown in Fig. 1. Hence the radii of the calculation lines 
are assigned the sign of the x coordinate, i.e., r is 
taken positive for compressors and negative for turbines. 
All values of angles are taken between -1/2 and +n/2. 


The coordinates of the first streamline approximation may 
be obtained by a linear interpolation of the calculation 
lines in cylindrical coordinates and then converting to 
Cartesian coordinates. Geometric parameters to be 
required in the analysis are calculated from the stream- 
line coordinates. The relative flow angle, 6, the 
distance along the streamlines, s, and the relative 
streamline curvature, K, are calculated at each calcu- 
lation point. The quantity (dA/ds)/A, where A is the 
area between and normal to the streamlines, is calculated 
for all calculation points except those on the blade 
surfaces. Required derivatives may be obtained using 
standard computer software for curvefitting and differen- 
tiation or other suitable methods, It has been found 
that good results may be obtained using a quadradic 
piecewise curvefitting algorithm to calculate the relative 
flow angle and thence the relative streamline curvature 
by dé/ds. 


The pressure, temperature and other fields of the rotor 
are functions of time relative to an absolute reference 
frame. Since they are not time-dependent, but are 
steady state relative to the rotor, the analysis is 
performed relative to the rotor with the exception of 
the energy and momentum analysis which must be related 
to the absolute frame, since the rotating reference 
frame is non-Newtonian. The equations are therefore 
considered applicable instantaneously, that is, variations 
with respect to time in the absolute reference frame are 
not considered. It should be noted that the static 
pressure and static temperature are properties of the 
fluid and are thus the same in the rotating and absolute 
reference frames, Whereas the greater part of the 
analysis is performed relative to the rotor, quantities 


which may be either absolute or relative, such as pressure, 


temperature, velocity and Mach number, will only be 
subscripted A when absolute. Total conditions will be 
denoted by capitals and static conditions by lower case 
letters. 


The following aerodynamic quantities are taken as 
specified: inlet relative total temperature, T,; inlet 
relative total pressure, Pj; mass flowrate, m; angular 
velocity of the rotor, w; angular acceleration of the 
rotor, w; specific gas constant, R; ratio of specific 
heats, y; and specific heat at constant pressure, cp: 
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; 
; 


The distribution of relative total pressure and relative 
total temperature is determined by relating downstream 
values to values at the initial or previous calculation 
lines. This is done by applying the Euler turbine 
equations as follows. The equation of the conservation 


of energy between two points along a given absolute 
streamline is 


“ ‘< 
Spt. +3 ie Rot = jdt, [2] 


where E is the energy transferred into (positive) the 
fluid per unit mass. The Euler turbine equation, which 


is a form of the equation of conservation of impulse 
and momentum, is 


E = (UV, sina) 547 Coes May sina); , (3] 
Noting that V, sina = U + V sing, and [4] 
v2 = v2 + u? + 2uv sins, [5] 


where a and 8 are the absolute and relative flow angles 
with respect to the radial direction, respectively, It 
can be shown by combining the THA and the Euler 


turbine equation and eliminating Ve and V, Sina that 


2 
(ce +42 - of silaee gevtouns? 
P 2 2 Ji Pp 2 Zz / j+1(6] 


which shows that the quantity 


2 2 
Vv U 
t+ "=, 
fp ee 


is constant along a given streamline. Further, noting 
that 


y2 
T=t+— and U = or {7], [8] 
2c 
P 
and that the initial relative total temperature is 


specified, the relative total temperature at any radius 
may be determined by 


ez ee 
Mecleny (tye wang |). ee ae u [9] 
2c, 
p 2c) 
40 


Whereas the flow is assumed to be ideal, adiabatic, and 
the relative total temperature is a function of a radius 
only, the relative total pressure is also a function of 
radius only and it may be determined by 


ne = 
Es PiledWcisdva kactaoa ue 
j-1 


The static pressure and static temperature distributions 
along calculation lines are obtainable by means of radial 
equilibrium analysis. This consists of simultaneously 
solving, along a calculation line, the Euler equation 

of motion and the continuity equation having relative 
total pressure, relative total temperature, and geometric 
parameters as known quantities. The Euler equation of 
motion relating the static pressure and acceleration 
fields applied to a particle of fluid at a given instant 
is 


SSSSe Vp. Salle « [11] 
9 


The acceleration of particle of fluid under consideration, 
a,, may be formulated mathematically in several ways. 
This analysis is based on that formulation of a, given 

by equation [1]. The magnitudes and direction of the 
acceleration components of a, are as follows: The 
acceleration component a, is further_seperable into two 
physically identifiable Components, a*, and @s,. 

ax, is the component which is directed towards the 
center of curvature of the relative streamline and is 
normal to the relative velocity. The equation for ay 
is obtainable by applying Euler's equation of motion 
to the motion of a particle in a plane expressing the 


gradient operator in cylindrical natural coordinate form 
and is 


n 


[12] 


as. is in the direction of the relative velocity. The 
eqUation for 4,_ is obtainable from the analysis of one- 
dimensional isentropic flow through a channel of 
varying flow area and is 


“ ee 
Be Uke tT, = m2 7a \A as Us [13] 


41 


The two foregoing equations are in terms of relative 
parameters, although the analysis for obtaining them is 
performed as if the reference frame were static. This 
is consistant with the concept of the acceleration 
gomponent a@%. Noting that = wi, V = Va,, fF = ri,, 

& = ga, and the definitions of the unit vectors as 
per Fig. 1, the equations for the Coriolis, centripedal 
and tangential acceleration components are 


For = 2uxV = 2uVun, [14] 
Gop = ox(oxF) = -wrt,, and [15] 
By = txF = iri, [16] 


respectively, The centripetal and tangential accelera- 
tion may each be separated into normal and streamwise 
components as follows: 


acpn = w*r sing un paz) 
= é 2 - 

acps -w r cosp U, [18] 
arn = wr cosé u, [19] 
@r, = or sing u, [20] 


Adding acceleration components having like unit vectors 


yields 
for Us, an, + Bcps + ars = as, and [21] 
for u,, Aan * Acpy * Acor * Bry [22] 


Substituting the foregoing into the Eurler's equation 
of motion and equating terms of like unit vectors yields 


ap > wesed 
x ua, pa, [23] 
AP g Seno a 
eal u, pal. [24] 
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The differential of static pressure along a calculation 
line between two consecutive streamlines is 


ee) pk 
dp = oe ds + mn dn. [25] 


The differentials ds and dn, in terms of dx and dy are 
ds = dx cosé + dy siné, and [26] 
dn = -dx siné + dy cosé, [27] 


The static pressure at an intersection point of a stream- 
line and a given calculation line is obtainable in terms 
of the static pressure at the intersection point of the 
previous streamline and the same calculation line by 
substituting the appropriate terms of the foregoing 
equations in 


Pi+t,j = Pig * - [28] 


In order to begin computing the static pressure distribu- 
tion along a calculation line, the static pressure at 
the one calculation point on the same line must be known. 
Whereas static pressures are not specified, a static 
pressure or parameter enabling the calculation of a static 
pressure must be assumed, It is best to assume the 
relative Mach number at the initial calculation point on 
a calculation line and calculate the SEAR pressure by 
y=) 


p=P (1+ M2) [29] 


It has been found desirable to choose the central stream- 
line, rather than the streamline indexed 1, as the 
initial streamline on which the Mach number is assumed 
and integrate outward toward the blade surfaces. The 
static pOnnereeuze, relative velocity and density at that 


point may be calculated by 
yl 
= TVG a) s [30] 
V=M YyRt, and [31] 
e = p/Rt. [32] 
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The static pressure differential and the static pressure 
on the next point on the calculation line may now be 
calculated using equations [12] to [29] and 


fi ee 


The relative Mach number at the next point is now known 
and the procedure may be repeated from point to point 
until the static pressure, static temperature and 
relative velocity distributions along the calculation 
line are obtained. Using the continuity equation 


m = SpVdA, [34] 


the total mass flowrate is obtained by integration along 
the calculation line. If the calculated mass flowrate is 
not equal to the specified mass flowrate, m » within 
an acceptable tolerance, a new estimate of Re relative 
Mach number is assumed, noting that 


P 
ee 1 [35] 
=i tae? 
where £(M) = M(1 + mM’) (7-1) 
TIRE 
oe £ yey" £0 15 oD — 


The values subscripted "OLD" are those based on the 
previous estimate of Mach number. This procedure is 
repeated until the calculated mass flowrate lies within 
an acceptable tolerance. This procedure yields the 
simultaneous solution of the continuity equation and 
the Euler equation, If the positions of the streamlines 
used are not a good approximation, the mass flowrate 
between adjacent streamlines will not be equal. New 
positions of the calculation points are computed such 
that the mass flowrate between streamlines is equal. 
The entire procedure is repeated successively until the 
repeatability of the parameters at each calculation 
point lies within an acceptable tolerance. The parameters 
selected for the repeatability check may be fundamental 
parameters, such as streamline positions, or required 
parameters, such as static pressure, static temperature 
and relative velocity. The root mean square of the 
relative velocity change between successive iterations 
is found to be most practical. This quantity may also 
be used for a convergence check between successive even 
iterations. Repeatability of parameters signifies an 
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approximation to the solution of the governing aerodynamic 
and geometric equations. The iterative process of deter- 
mining the streamline positions may diverge if the 
previous approximation of the positions is a poor one or 
if the change between successive iterations is too large. 
The use of "damping factors," which limit the changes of 
the streamline positions, curvatures, and the quantities 
dAds/A, prove useful in attaining convergence. When 
using compas factors, a minimum number of iterations, 

4 to 6, should be allowed before testing for repeatability 
to ascertain that repeatability is due to the attainment 
of a solution and not to excessive damping. It is further 
desirable to increase mass flowrate and angular velocity 
of the rotor in steps until desired values are attained. 


The accuracy of a numerical solution may be assessed in 
several ways. For flows in static channels, a numerical 
solution may be compared with a readily available 
analytical solution, such as a solution for the flow in 
a channel bounded by logarithmic spirals. Whereas 
analytical solutions for flows in rotating channels are 
not readily available, the accuracy of a numerical 
solution may be assessed by comparing the power input or 
output calculated by three different methods: by energy, 
by momentum and by integrating the static pressure on 
the blade surfaces. Whereas the law of conservation 

of energy and the law of conservation of impulse and 
momentum are applied using data along inlet and outlet 
calculation lines and the static pressure integration 
procedure is applied using data along bounding stream- 
lines, good agreement between the three values of the 
power tends to indicate a good solution, although this 
does not constitute an absolute proof, The calculation 
of the torque about the centerline of rotation by two 
different methods may also be used for assessing the 
accuracy of solutions for both static and rotating 
channels. 


Whereas the anlysis assumes adiabatic inviscous ideal 
flow--i.e., no losses, boundary layer, secondary flows-- 
the solution is one for ideal primary flow only. Under 
most circumstances the primary flow will be perturbed 

by secondary effects; however, primary flow character- 
istics will predominate. As the magnitude of the 
perturbations increases, secondary effects may eventually 
predominate, e.g., flow separation. However, the primary 
flow does indicate basic characteristics of the flow in 
turbomachines and may serve as a basis for analysis of 
non-ideal flow. Probable boundary layer migration and 
flow separation may be deduced from the static pressure 
and relative velocities fields of the ideal flow. It 
should be noted that total pressure losses in real flow 
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will affect the nature of the flow. Methods of including 
total pressure losses in the analysis are currently being 
investigated by the author. It should also be noted that 
boundary conditions have been implicitly determined by 
the algorithm for determination of relative streamline 
curvatures, This algorithm yields good results as 
indicated by the power criteria outlined. A detailed 
analysis of boundary conditions is beyond the scope of 
this paper. Such an analysis must consider the flow 
outside the rotor. Attempts by the author at forcing 
compatability between the flows inside and outside of 
the rotor by forcing uniform inlet, or both inlet and 
outlet, relative velocities have been found to yield 
poor results as indicated by the power criteria outlined. 
It is the opinion of the author that the interaction 
between the flows inside and outside of the rotor is 

one of the justifications for the assumption of circum- 
ferentially constant inlet relative total temperature 
and pressure. Circumferential variations in inlet 
relative total temperature and pressure may be based on 
circumferentially constant inlet absolute total tempera- 
ture and pressure, Some references on experimental work 
which give some insight into real flows are included in 
the reference section. 


Blade surface relative velocity distributions and 
relative streamlines calculated for a rotor having 45 
degree logarithmic spiral blading are shown in Fig. 2. 
This and other such data obtained from a program based 
on the method of analysis outlined in this paper is 
useful in the discussion of the mechanics of turbo- 
machinery from both the theoretical and practical point 
of view. The analysis, although based on equations 
usually studied in elementary courses on machines and 
on fluids, is similar to the type of analysis which 
might be performed by one in industry and may help 
bridge the gap between textbook and "real-world" problems. 
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NOTE: A computer program for the above analysis using 
a CDC computer may be obtained by writing to the author. 
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